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We present the inconsistencies which arise in quenched QCD in A J = 1/2 non-leptonic kaon decays using chiral 
perturbation theory (xPT) to one loop. In particular we discuss how the lack of unitarity of the quenched theory 
invalidates the usual methods for the extraction of matrix elements from correlation functions. 
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1. INTRODUCTION 

While different techniques have been devel- 
oped in order to extract A/ = 3/2 K tttt 
matrix elements from the lattice and fi- 

nite volume corrections can be controlled even 
in the quenched approximation (at least to one 
loop in xPT) (|,|^,|,|), this is not true for the 
A/ = 1/2 channel. In this case, in addition to 
the presence of large final state interactions (FSI), 
there are manifest inconsistencies due to quench- 
ing. The reason why the A/ = 3/2 case is dif- 
ferent will be explained in sec. |[ while the ex- 
plicit calculation of the scalar form factor (which 
shares the same quenching artifacts as K tttt, 
AI = 1/2) is in ref. [Q. Some of these problems 
in the xPT framework were first discovered in 
rcfs. |^,||j9) , our aim is to exhibit the full set of the 
quenching effects and to understand their origin. 
By computing matrix elements involving isoscalar 
TTTT scattering contributions (e.g. K tttt or the 
scalar form factor) using quenched xPT (QxPT), 
in finite and infinite volume, with the insertion of 
different AI =1/2 weak operators we observe 
that: 

1) The strong phase shift is no longer universal 
since it depends on the choice of the weak opera- 



tor. 

2) The Liischer quantization condition and the 
Lellouch-Liischer (LL) finite-volume correction 
factor for matrix elements are no longer valid. 

3) The standard procedure for extracting 
matrix elements from time-independent ratios 
of finite-volume euclidean correlation functions 
(CF) fails, since the time dependence of the CF 
depends on the choice of the two-pion sink. 

4) The 77' double-pole produces new polynomial 
terms in t (time) and L (box size) which cause 
practical problems in the extraction of matrix el- 
ements. 

2. FULL (UNQUENCHED) QCD 

The unquenched K — > tttt matrix element can 
be extracted by studying the time behaviour of 
a suitable four-point euclidean CF in finite vol- 
umes: 
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J2{0\TT^{0)7rt^.mW){W\HwmK) X 

x(/^|/^l(0)|0)e-™^l*^l-'^l*l , (1) 

where (pqit) = J cfix(l){x,t)e'^'^'^ , \W) is a two-pion 
state with total energy W, H,^ = {q : \q\ fixed} 
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the matrix element is extracted |^) and finally 
the last two terms represent the chiral expansion 
of the LL finite volume factor to this order. The 
validity of eq. (||) relies on the unitarity of the 
theory. 



Figure 1. The one-loop diagram which contains 
power corrections in the volume. 



and the large time limit has been taken while 
isolating the contribution with time behaviour 
^-w\t\_ jji order to extract the matrix element 
{W\Hw{Q)\K) one has to divide by the fac- 
tors corresponding to the kaon source and the 
two-pion sink extracted from the corresponding 
CF H. The validity of eq. (|l]) and Watson's the- 
orem rely on unitarity, and it is the absence of 
unitarity in the quenched theory which is respon- 
sible for most of the difficulties listed at the end 
of sec. ^. The calculation in infinite volume is 
straightforward and we will not discuss it further, 
except to point out that of the one-loop diagrams, 
it is only that in fig. which can have physical 
intermediate states and hence an imaginary part 
(and 1/L" corrections in finite- volumes). 

2.1. Finite volume 

Evaluating eq. (Q) to one loop in xPT on a fi- 
nite volume (FV) one obtains, up to exponentially 
small terms in the volume. 



^-mK\tK\ p~2E\t\ 
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where Aoo is the corresponding infinite volume re- 
sult, T parameterizes FV corrections, AW is the 



energy shift (cx l/L^), 



and z(i) 



are defined for example in [§,1). The first two 
terms in T shift the two-pion energy in the time 
exponential and in the argument of Aoo respec- 
tively, the third term depends on the sink used 
to annihilate the two pions (it is cancelled when 



3. QUENCHED QCD 

We now consider the determination of the ma- 
trix elements in quenched QCD, restricting our 
attention to the diagram in fig. (|l]), which is the 
one leading to anomalous effects. In particular it 
receives two unphysical contributions: (i) ghosts 
which cancel internal quark loops, (ii) 77' with 
double-pole insertions (^ ) . They lead to a number 
of inconsistencies: 

(i) The ghosts are a manifestation of the lack 
of unitarity. Phase shifts, energy shifts and fi- 
nite volume corrections are no longer universal, 
the factor corresponding to the sink cannot be 
removed and the extraction procedure in eq. (|l|) 
cannot be applied. 

(ii) rj' double-poles lead to a more severe and 
unphysical infrared structure. In infinite volume 
they appear as divergences in the chiral limit (the 
so-called "quenched chiral logs", see for exam- 
ple ||]) and in the singular limit at threshold 
|0,|lO I . In finite volumes they produce terms with 
an anomalous behaviour in time and volume (see 
below and [0|,|l). 

3.1. The scalar form factor 

To study these effects consider the time-like 
scalar form factor of a scalar and isoscalar op- 
erator S, which shows all the effects of if ^ tttt, 
AI — 1/2 but has a simpler chiral structure. The 
corresponding CF at one loop in QxPT can be 
written as follows: 

„-2E\t\ 

(0|4(i)7r:,-(i)^(0)|0) - ^^^^[1 + T] , (3) 
ait + asi^ + asi^ 60 + bit + b2t^ 



T = 

{ELf 
coz{Q) ciz{l[ 
' [ELf EL 



{ELf 

+ C2z{2){EL) + C3z{3){ELf 



^ Both these contributions do not appear to one loop in 
the A7 = 3/2 case because of total isospin conservation. 
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Notice that the term dAac/dW is not present 
because in this case Aoo is a constant. The ai 
term would be the quenched energy shift, but it 
now depends on the choice of operator S and re- 
ceives contributions from rj' double-poles. The 
remaining a^'s come from the rj' double-pole in- 
sertions giving an anomalous behaviour in time. 
The terms proportional to the &i's depend on the 
sink and are different from the corresponding ones 
in full QCD. A consequence of the lack of unitar- 
ity is that they cannot be eliminated using the 
procedure described in sec. ^. The terms propor- 
tional to the Ci's correspond to the finite- volume 
corrections. They are no longer universal and the 
last two-terms, which come from the rj' double 
pole, appear prima-facie to diverge in the large- 
volume limit. We now discuss this point further. 

4. VOLUME DEPENDENCE 

The apparent volume dependence in eq. (^) is 
rather strange. By evaluating the matrix element 
in infinite volume and euclidean metric no diver- 
gences appear away from the chiral limit, and one 
would hence expect that the large volume limit 
of the finite volume matrix element should ex- 
ist. In spite of appearances this could be the 
case, since at fixed physical momenta q (so that 
the exitation level n grows as the volume is in- 
creased), the z(j)'s scale with the volume, for in- 
stance z(0) — —V (X (Q). By dimensional anal- 
ysis and numerical evaluation, it is possible to ver- 
ify that the z(z)'s decrease with L fast enough to 
compensate the powers of L appearing in eq. (||) . 
However we do not yet have a proof of this since it 
is not possible to perform the large volume limit 
at fixed physical energy without the quantization 
condition. 

5. CONCLUSIONS 

We have shown that quenching artifacts in 
K TTTT, A/ = 1/2 matrix elements are very 
important and might invalidate the usual proce- 
dure for extracting physical matrix elements from 
lattice CF. Problems related to the non-unitarity 
of the quenched theory also affect K tttt ma- 

"where u = J^n^ ""^^^ = {n e : |27rn/L| = 



trix elements with AI — 3/2, but they do not 
appear at one loop in xPT and are therefore less 
severe. 

The problems we have raised do not necessarily 
have a solution within the quenched approxima- 
tion, which, it must be remembered, is an intrin- 
sically unphysical theory. However, it is worth 
exploring whether there is some practical way for- 
ward, which is not obviously inconsistent. 

Quenched QCD is perturbatively reproduced in 
the xPT framework by enlarging the set of parti- 
cle states to contain new unphysical mesons. This 
extension can be realized in two different, but per- 
turbatively equivalent, ways either with the su- 
persymmetric enlargement of the chiral group |^] 
or with the replica method We are cur- 

rently investigating properties of the enlarged set 
of states and the behaviour of weak matrix el- 
ements between those states, including the un- 
physical ones 0. In particular, we are exploring 
whether it might be possible that some of those 
universal relations of full QCD mentioned in sec. 
can be generalized to the quenched approxima- 
tion, once the extended set of external states is 
considered. 

REFERENCES 

1. SPQCDR Collaboration, Ph.Boucaud et al., 
Nucl. Phys. Pro c. Suppl. 106 (200 2), 323 

2. C.-J.D. Lin et ai, |hep-lat/0208007 

3. L. Lellouch and M. Liischcr, Commun. Math. 
Phys. 219 (2001), 31 

4. C.-J.D. Lin, G. MartineUi, C.T. Sachrajda, 
M. Testa, Nucl. Phys. B619 (2001), 467 and 
Nucl. Phys. Proc.SuppL 109 (2002), 218 

5. C. Bernard and M. Goltcrman, Phys. Rev. 
D46 (1992), 853 

6. M.F.L. Golterman and K.-C. Leung, Phys. 
Rev. D56 (1997), 2950 

7. C.-J.D. Lin et ai, in preparation 

8. C. Bernard and M. Golterman, Phys. Rev. 
D53 (1996), 476 

9. M.F.L. Golterman and E. Pallante, Nucl. 
Phys. Proc. Suppl. 83 (2000), 250 

10. G. Colangelo and E. Pallante, Nucl. Phys. 
B520 (1998), 433 

11. P.H. Damgaard and K. Splittorff, Phys. Rev. 



4 

D62 (2000), 054509 



